Anomalous tricritical behaviour in the coil-globule transition of a single polymer chain 
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We investigate a model of self-avoiding walk exhibiting a first-order coil-globule transition. This first-order nature, unravelled 
through the coexistence of distinct coil and globule populations, has observable consequences on the scaling properties. A 
• thorough analysis of the size dependence of the mean radius of gyration evidences a breakdown of the plain tricritical scaling 
behaviour. In some regimes, anomalous exponents are observed in the transition region and logarithmic corrections arise along 
the coexistence curve. 
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Experimental observations ||] 13] as well as theoreti- 
cal studies for a specific model [||| Q] suggest that the 
coil-globule transition of an isolated self-avoiding walk of 
finite size N might exhibit first-order features. A strong 
evidence is the coexistence in some range of tempera- 
tures of distinct coil and globule populations. Its theo- 
retical signature is a bimodal shape of the order param- 
eter distribution pj^\t), where t is some conformational 
parameter of the single chain, such that the transition 
is achieved through an exchange of weight between the 
two peaks as temperature varies. We here investigate the 
observable consequences of this first-order feature on the 
scaling behaviour of the chain in the transition region. 
Although one recovers a second-order transition in the 
infinite-size limit, one may wonder whether the standard 
tricritical picture of the 0-point |^ still applies. 

In previous studies |^ Q, scaling arguments supple- 
mented with Monte Carlo simulations of a self-avoiding 
chain on a cubic lattice provided an exact analytical ex- 
pression for the entropic contribution Pj^~^\t) to the 
distribution Pj(f (t) of the order parameter t, defined as 
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where r is the radius of gyration of the conformation and 
p = N/r'^ its mean density (here d = 3 and = 3/5 is 
the Flory exponent). We then modelled the energy of the 
isolated chain a,s U = —NJt. Scaling analysis detailed in 
previous studies H M led to introduce rescaled variables: 
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f = tN'^ with 4>=\-- and r = 1 - - (3) 
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where 6* is a rough estimate of the transition temperature. 
We obtained the following expression for the equilibrium 
distribution: 



Pjv(f,f) = 



h{f, t) e 



_yl'JV-5r(l-l/")4-? 



UN,t) 



(4) 



where 
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is the scale-invariant contribution. The factor Xc(iV, f) 
ensures that PN{f,i)di = 1. Let us emphasize the 
key role of the factor t'^: for c < —1, this factor, hence 
the function /i(f, t), are not integrable in i = 0. It forbids 
to take the infinite-size limit in Pi^[t,t) and to use the 
infinite-size system as a reference system. This factor if 
outweights the coil region and it will induce an observable 
breakdown of the scale invariance. More generally, we 
introduce for any real z: 
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An arbitrary moment simply writes: 
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Due to the factor , the scaling behaviour of Iz {N, f) 
strongly depends on the sign of z + 1. For z + 1 > 0, it 
writes: 
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In the case where z + 1 > 0, Tz{N,f) is thus asymp- 
totically scale- invariant, depending only on the rescaled 
variable f = N'^t but no more on N. 

For z + 1 < 0, two regimes occur on each side of a 
borderline Tcoex{N, z), as shown on Figure 1: 
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Here arises the above-mentionned scale- invariance break- 
down. It originates in the bimodal shape oft^^'^pNif, t): 
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when f decreases, a "coil" peak is quite abruptly re- 
placed by a "globule" peak after a coexistence regime 
precisely located around the curve Tcoex{N,z). Writing 
that the two peaks of t^'~'^Piy{f,t) coexist with equal 
weights (which requires z + 1 < 0) yields the following 
implicit equation for the borderline Tcoex{N, z): 
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where is some numerical constant. Tcoex{N,z) could 
have been obtained straightforwardly by writing the bal- 
ance between the two expressions of Iz given in Eq.(||). 
For large A'", Tcoex{N, z) oc (InTV)*^ for any value z < —1. 
The borderline ■Tcoex{N,z) thus separates the domain in 
(f, lniV)-space where the globule peak is overwhelming, 
leading to a scale- invariant expression for Iz , and the do- 
main where the coil peak is overwhelming. In the latter 
domain, the behaviour of Iz is ruled by the coil peak, 
being controlled only by the size N (independent of f at 
fixed N). Note that we are here reasoning on the peaks 
of i^~'^PN{T,t) so that the two peaks cannot be exactly 
identified with the coil and the globule phases, unless 
z = c. 

In the special instance where z = c, fcoex (N^ c) is then 
the "true" coexistence line where the coil and globule 
populations are equally weighted; it writes more explic- 
itly (denoting 6 — 9{oo)): 
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Let us note that the value n = 2 is consistent with the 
predicted value 4> = 1/2 for the crossover exponent ||] 
[|6[, since here ~ 1 — 1/n. 

Observable quantities are related to the moments 
of the distribution PN{T,t). Recalling that r — 
j^i/d^i/d~u ^ j^u^ii/d-^ where: 
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we introduce a generalized mean radius of gyration as: 

= A^i/d < p >-^^^ TV"" < r >-^^ (13) 

for any nonzero real a. Scaling behaviour of the mo- 
ment < > follows from the behaviour of Ic and 
la+c (see Eq.(^). Of special interest are the values 
a = 1, as i?i is related to the mean value < t >, and 
a — — 2(z^ — 1/d), leading to the standard mean radius of 
gyration Ro =< >^/^, that can be measured through 
scattering experiments. For fixed r < 0, one recovers the 
foreseeable globule scaling law Ra oc N^^'^ (here d = 3) 
whereas for fixed r > 0, one recovers the coil scaling law 
Ra oc N'^ , whatever a is. In between, in the transition 



region |r| — > 0, f finite, various scaling regimes are ob- 
served, depending on the sign of 1 -|- c and 1 + c + a. The 
different cases are sketched on Figure 2. 

If c > — 1 and a + c> —1 (case 1), the ratio la+c/Ic is 
scale-invariant; it follows that Ra exhibits the standard 
tricritical behaviour: 
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For negative and large f — N'^t, one recovers the usual 
scaling behaviour Ra oc N^^'^; this crossover reflects in 
the hyperscaling relation: 
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whereas for positive and large f = N'^'t, one recovers the 
usual scaling behaviour Ra oc N"; this crossover reflects 
in the hyperscaling relation: 
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As soon asc<— lorQ!-|-c<— 1 (cases 2 to 5), scale- 
invariance breakdown and anomalous behaviour are ob- 
served. The scaling law for Ra still writes: 
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But the exponent I'eff and the scaling function / 
change when passing across the curves Tcoex [N, c) and 
Tcoex{N, a+c). It is clear that (/) = 1 — 1/nis the crossover 
exponent. In any cases (2 to 5), both la+c and Ic are 
scale-invariant in the leftmost region. For f < with |f| 
enough large, any observable Ra thus exhibits the scaling 
behaviour yet encountered in case 1: 
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For |f I ^ 1, one recovers the usual scaling behaviour 
Ra oc A^^/'', for any a, according to the hyperscaling 
relation given in Eq. (p^ . 

In cases 2 to 5, the key point is the emergence of an 
intermediate region where either la+a either Xc, or both 
behave as f-independent powers of N . This modifies the 
exponent vg into an anomalous exponent: 
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observed for a < 0. The axis f = plays no particular 
role. In cases 2 and 3, only one curve exists, respec- 
tively Tcoex{N,a -I- c) and Tcoex{N,c). On its right side 
(i.e. above it), Ra behaves respectively as N'^f> f2{T) and 
N'^" fz ij) where /2 and /s are complicated (but indepen- 
dent of N) functions of f . Both reduce to a power law 
for f > 1: 
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Using definitions ( |20D and (|2l|), we obtain hyperscaling 
relations: 
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ensuring that the expected scaling behaviour oc N'^ 
is recovered as soon as f ^ 1. Note that the mean or- 
der parameter < t > and the associated mean radius of 
gyration i?i belong to case 3 (c = —1.13 and a = 1). 

In the last cases (4 and 5), the anomalous exponent 
is observed in the intermediate but unbounded region 
located between Tcoex {N, c) and Tcoex (N, a + c). It is still 
equal to ly'gia) for a > and i^g{ct) for a < 0. In these 
cases, the exponent v is observed in the rightmost region, 
even if f is not large with respect to 1 (even negative). 
Note that the mean radius of gyration < >i/2 belongs 
to case 4 (c = — 1.13 and a + c = —1.66). 

The remarkable features unravelled here are: 

(i) the exchange of weight between the coil phase and 
the globule phase is far sharper with respect to the vari- 
ation of f when it occurs through the coexistence of two 
peaks than when it occurs through the shift of a simple 
peak. Accordingly, a sharp crossover is observed in the 
behaviour of if z -I- 1 < 0, located in a narrow stripe 
of width (logA^)""^ around Tcoex{N, z), whereas a slow 
change with a scale-invariant behaviour of is observed 
for z -I- 1 > 0. 

(ii) due to the A^-dependence of the curves Tcoex{N, z) 
for z -I- 1 < 0, a size-controlled crossover occurs at the 
passage across the curves Tcoex {N,c) (if c < —1) or 
TcoexiN,a + c) (if Q 4- c < -1). 

(iii) when either c, either a -I- c or both are lower than 
— 1 (cases 2 to 5 of Figure 2), an intermediate region 
arises where an anomalous exponent (neither v nor vg) 
is observed. Usual tricritical scaling is observed only for 
c > — 1 and a + c > —1 (case 1 of Figure 2). 

(iv) for c < — 1 (cases 3 to 5 of Figure 2), on the co- 
existence line TcoexiN,c), the scaling behaviour of Ra 
exhibits multiplicative logarithmic corrections with the 
same exponent whatever a is: 



What we here suggest is that the scaling behaviour 
around the 0-point might be more complex than the 
standard tricritical behaviour Q The peculiarity 

of the scaling behaviour associated with the first-order 
features of the transition is the occurence of two "theta 
regimes" . A first one is observed when the globule peak 
is overwhelming; it thus involves the exponent vg arising 
from the scale invariance of the globule phase and the 
scaling function is a power law. A second, anomalous 
one involves an observable-dependent exponent i/g (a) (if 
c < — 1 and a > 0) or Vgia) (if a -I- c < —1 and a < 0) 
where a refers to the observable Ra', this anomalous ex- 
ponent arises from a non trivial balance between the scale 
invariant globule phase and the size-controlled coil phase, 
which is required to match their incompatible scaling be- 
haviours. Only the detailed analysis of equilibrium distri- 
butions hke Pj^\t) gives reliable predictions about the 
actual coil-globule transition experienced by a chain of 
finite size. 

A related result is the irrelevance of standard finite- 
size scaling approaches to describe the size dependence of 
the thermal coil-globule transition. Indeed, the standard 
finite-size analysis used for first-order transitions Q here 
failed as the two states merge in the infinite-size limit. 
On the other hand, an analysis based on the knowledge 
of the infinite-size second-order transition cannot ac- 
count for the first-order features here observed in finite 
size. In a word, standard finite-size scaling approaches 
are designed to predict the rounding of the transition fea- 
tures in finite size but cannot introduce back the change 
of nature of the transition observed here. A novel finite- 
size scaling analysis should thus be designed to interpret 
experimental or numerical data [pil. 
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Captions 

Figure 1: Scaling behaviour of Zz for z -|- 1 < (here 
z = —1.66, corresponding to the integral la+c involved in 
the mean radius of gyration < >^/^). The two scaling 
regimes are separated by the curve TcoexiN, z), depending 
on z but of similar shape for any value z < —1. The curve 
is restricted to the domain f < ig w —3.2 (independent of 
N) where a well-identified globule state (a globule peak) 
exists. For c < — 1, fcoexiN. z = c) is the coexistence line 
of the then first-order coil-globule transition. 

Figure 2: Critical exponents of Ra(N, f) for the differ- 
ent scaling regimes determined by the signs of 1 -I- c and 
l+a+c. The bold line (cases 3,4 and 5) is the coexistence 

curve Tcoex 

{N,c) (defined for c-M < 0, here c = -1.13); 
the thin line (cases 2, 4 and 5) is the curve fcoex{N, a + c) 
(defined for 1 + a + c < 0, here a + c = —1.66 which cor- 
responds to the mean radius of gyration < >^/-^). 
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Case 1: c < — 1 , q; + c < — 1 





